On the structural properties of small-world networks with finite range of shortcut links 
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We explore a new variant of Small- World Networks (SWNs), in which an additional parameter 
(r) sets the length scale over which shortcuts are uniformly distributed. When r = we have 
an ordered network, whereas r = 1 corresponds to the original SWN model. These short-range 
SWNs have a similar degree distribution and scaling properties as the original SWN model. We 
observe the small-world phenomenon for r<l indicating that global shortcuts are not necessary for 
the small- world effect. For short-range SWNs, the average path length changes nonmonotonically 
with system size, whereas for the original SWN model it increases monotonically. We propose an 
expression for the average path length for short-range SWNs based on numerical simulations and 
analytical approximations. 



I. INTRODUCTION 

The structural and dynamic properties of small-world 
networks (SWNs) [l[ have been studied intensively with 
applications to diverse fields such as sociology, biological 
networks, neural networks and so on [34Tl|. In the orig- 
inal model, the SWN is constructed by adding shortcuts 
on top of a connected lattice ring (ordered lattice). The 
added shortcut connects a pair of nodes picked randomly. 
This represents the case such that each node in the net- 
work has equal probability to make random connections 
to all others regardless of their separation for the ordered 
lattice (also denoted as Euclidean distance). Under this 
construction, the small- world effect appears in SWN with 
the mean separation between any pairs of nodes scaling 
logarithmically with the network size. 

Although there has been intensive work on the struc- 
tural properties of SWNs, the relationship between the 
maximum length of shortcut and the small-world effect 
has not yet been investigated, to the best of our knowl- 
edge. In the original SWN model, the maximum shortcut 
length is half of the network size. On the other hand, if 
we limit the maximum shortcut length to be some small 
value, e.g. corresponding to the case that the shortcut 
can only connect the nearest neighbors, the network will 
be identical to a regular lattice. These limits indicate 
that by tunning the maximum range of shortcut lengths, 
we can change the network form a standard SWN to a 
regular network. An important question that arises is: 
what are the constraints on the maximum shortcut length 
in order to give rise to the small-world effect? 

Besides the motivation stated above, there are other 
reasons to examine more closely the role of shortcut range 
in SWNs. In the standard SWN model, shortcuts are 
global with uniform length distribution, which assumes 
that the shortcut connection is completely random. How- 
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ever, this may not be an accurate representation for real- 
world networks. For example, in social networks it is 
reasonable to assume that the probability to have a con- 
nection is higher for two individuals living close by as 
compared to two individuals living far away. One simple 
way to address this issue is to modify the property of the 
shortcuts in SWN. In some previous studies, a variant of 
SWN was proposed with global shortcuts but power law 
distribution of the shortcut length [r3 - [h4T |. This repre- 
sents the case that the further apart the nodes are on 
the ordered lattice, the harder it is to have a shortcut 
between them. 

In this work, we study another variant of SWN such 
that the shortcut length has a uniform length distribution 
but with limited range (maximum length). We call this 
model the short-ranged small- world network (SRSWN). 
In SRSWNs, we introduce a new parameter r which rep- 
resents the maximum proportion of the network that a 
shortcut can reach. This corresponds to the case that one 
has equal probability to build random connection only 
within a community. Moreover, as the r value varies 
between and 1, the SRSWN converts from a regular 
lattice to a standard SWN. The investigation of SRSWN 
will also provide insight into how the small-world effect 
depends on the shortcut length. In the following, we will 
explore the structural properties of SRSWN with both 
numerical and analytical approaches. 



II. MODEL AND RESULTS 

The construction of SRSWN is similar to that of 
SWNs. We first start with a regular connected lattice 
with L — 2N nodes. Shortcuts are added with probabil- 
ity p, with totally x = Lp added. Two types of distance 
are introduced: the "Euclidean distance" is defined as the 
distance of two nodes along the backbone and the "min- 
imal distance" denotes the shortest distance after taking 
shortcuts. An outgoing shortcut will randomly connect 
a node within Euclidean distance rN, i.e. a node ran- 
domly chosen from rL nearest neighbors. Note that this 
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FIG. 1. plots of numerical results of £ vs L with p = 0.004 
and various r values. As shown, I changes nonmonotocially 
in SRSWN while it changes monotocially in SWN. When L is 
large enough, I varies logarithmically with L, indicating the 
appearance of small-world property. 



construction yields the same degree distribution as that 
of SWN (see Appendix) . In this work, we concentrate on 
the case of small p(p < 0.01) values. For each realization 
of network, we pick a random initial position and apply 
the breadth-first search algorithm to find the quantities 
analyzed. The simulation data presented are averaged 
over at least 10000 network realizations. 



A. Mean path length I and l(n) 

To check if SRSWN displays small- world phenomenon, 
we first check the relationship between mean path length 
£ (defined as the average distance along the shortest paths 
for all possible pairs of network nodes) and network size 
L by simulations. The plot of I vs L in SRSWN and 
SWN is shown in Fig[T] which displays some interest- 
ing features: (i) Initially £ for the SRSWN increases al- 
most linearly with L, which is similar to a regular net- 
work, (ii) For large values of L, £ for SRSWN has similar 
behavior as that of SWN, i.e. I scales logarithmically 
with L. (iii) In the intermediate region, J first reaches a 
maximum value and then decreases as L increases. This 
nonmonotonic behavior points to the surprising property 
that nodes in the network come closer to each other (on 
average) as the network size increases. This feature is 
not observed in either regular lattice or SWN, as £ in 
both cases increases monotonically with L. Furthermore, 
point (ii) noted above indicates that SRSWN does show 
small- world properties as long as the network size is large 
enough. 

To gain a further understanding of the numerical obser- 
vations noted above, we first study £(n), which is defined 
as the minimal path length between two nodes given that 
the Euclidean separation is n. Alternatively, £(n) can be 



considered as the number of time steps a walker needs 
to take (via the shortest path) from origin to a node at 
Euclidean distance n away. The relation between £(n) 
and £ is: 



n—N 



N ^ V 



(1) 



In the regular network, the walker can only walk 
through the arc of the circle such that £{n) — n. In SWN, 
because the shortcut range is global, taking a shortcut 
will completely randomize the position of the walker. As 
the shortest path usually contains at least one shortcut 
and the shortcut can land on any node with same proba- 
bility, it seems plausible that there will be a large number 
of nodes having the same probability to be visited in a 
given n steps. In previous work, it has been shown for 
large L that £(n) ~ n when n < 1/p and it quickly satu- 
rates to £ as n increases [H, HH . This indicated that the 
majority of the nodes can be reached within the same 
time steps on average, regardless of their original posi- 
tions. 

In SRSWN, however, £(n) will behave in a different 
manner from either the regular network or the SWN. 
First, a shortcut in SRSWN will bring the walker closer to 
nodes far away and save some time steps, indicating that 
£(n) will be less than n for most of the nodes. On another 
hand, due to the limited range of the shortcut, it cannot 
take the walker to an arbitrary node as in the SWN, 
indicating that £(n) will depend on n. Let's examine the 
case that r is small and focus on the behavior of £(n) for 
nodes away from origin (n = 0) and the extremal node 
(n = N). In particular, let us consider three nodes with 
Euclidean distance m, n 2 and 77,3 such that ri2 = rii + rN 
and n 3 = n 2 + rN. Because a shortcut has the maximum 
length rN, node 77,2 cannot be visited before node n\. 
Clearly we have that ~£(n 2 ) > £(m) and £(n 3 ) > T.(n%). 
Denote 5£\ — ?(n 2 ) — ~£{n\) and S£% = £(n 3 ) — £(n 2 ) as 
the extra time steps the walker takes in going from node 
m to n 2 and n 2 to 713, respectively. Because the network 
configuration for the interval [m, n 2 ] is similar to that for 
[n 2 , n 3 ], it is reasonable to assume that the walker will, on 
average, take the same number of time steps in traveling 
between these two regions, i.e. 52% = 5£ 2 . Given that 
the Euclidean separations for the pair of nodes n\ ,n 2 and 
n 2 ,ri3 is the same (rN), we expect a linear relationship 
between £ (n) and n for <C n <C N. 

The simulation results for £(n) shown in Fig [2] confirms 
the expectations from the rough argument given above. 
As displayed, £{n) has a part that increases linearly with 
n, which corresponds to the interval n G [rN, (1 — r)N]. 
Furthermore, we notice that for networks with different 
size but same parameter p and R = rL, £(n) increases 
with n at the same rate. This is because the time steps 
the walker takes to go through the interval [711,^2] and 
[77,2,713] only depends on the chance of encountering a 
shortcut (respectively p) and the range of the shortcut 
(respectively rN or R). As the result, the rate corre- 
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FIG. 2. In the main figure, we show numerically computed 
t(n) with different network size L but same R and p values. 
Apart from the nonlinear behavior near n ~ and n ~ N , 
£{n) increases linearly with n in the middle region, presum- 
ably corresponds to the interval n £ [rN, (1 — r)N]. The slop 
of the linear part in the three cases are the same, indicating 
that the linear increase rate of l(n) depends on R and p. In 
the inset, we show the linear relationship between £ and L 
with the same choice of p and R values. The points can be 
well fitted by the line y = 450 + 0.04a;. 



sponding to the linear increase of £(n) depends on p and 
R = rL only. 

In the limit that R <C L (equivalently r -C 1) such 
that the nonlinear behavior of £(n) near the beginning 
and the end can be ignored, we can approximate the 
linear relationship between £(n) and n for all n's from 1 
to N. Under such an approximation, t derived based on 
Eq.(pQ) will be a linear function of L for given R and p 
values. It is noteworthy that while this approximation is 
for R <C L, we found from simulations (as shown in the 
inset of Fig [2]) that the linear relationship holds even for 
R values close to L. Given such linear relationship and 
considering that the minimum value of L for given R is 
L = R, we can express the £ of SRWN with parameter 
L, p and r as 

£(L, p, r) = (L — R)K(R, p) + £ SWN (R, p) , (2) 

where K depends on the linear increase rate of £(n) that 
is a function of R = rL and p. £swn(R,p) denotes the 
average distance of SWN with size R and shortcut prob- 
ability p, since the network that corresponds to L = R is 
the standard SWN. 



B. Scaling properties 

Eq.© provides a scaling form of 1 in SRSWN. On the 
other hand, the scaling properties of SWN have been de- 
rived from renormalization group transformations [l7l . ll8j 



and this approach can also be generalized for SRSWN. 
The renormalization group transformation consists of 
conserving the number of the shortcuts while varying net- 
work size L and shortcut probability p. As the parameter 
r is the length scale of the shortcut range, £ in SRSWN 
will continue to be an homogeneous function with r as 
another variable: 



which gives another scaling form of £ as: 
?{L,p,r) = Lg(x,r), 



(3) 



(4) 



where g is a scaling function and x — Lp is the total 
number of the shortcuts. The scaling relation in Eq.Q 
has been confirmed by our simulations. 
Combining Eq.(j4]) and Eq.((2]), we have: 



-=g(x,r) = (l-r)K(R,p) + 
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(R,p) 



(5) 



Furthermore, based on the scaling property of SWN [17], 
UM, we have 



-SWN 



(R,p) 



R 



= f(Rp) = f(rx), 



(0) 



By combining Eq. ([5]) and Eq. © , we can express the func- 
tion g(x, r) as 



g(x, r) = (1 — r)K(R.p) + f(rx)r, 



(7) 



which indicates that K(R,p) is also a function of x and r. 
By keeping r fixed but choosing different p and L value 
as p' — p/X, L' — LA, we keep x and r the same. This 
leads to the following relationship 



K(rL,p) = K(x,r) = K(XrL, ^). 

A 



(8) 



The only way to have Eq.® hold is that K is a. function 
of rLp only. Eventually, we have the simplified scaling 
form of £ from Eq.© as: 



£ = (L-R)K(y)+Rf(y), 



(9) 



with two scaling functions K and /. The variables are 
y = rx = rLp and R — rL. 



C. Numerical analysis 

The functional form of f(y) in Eq.© has been stud- 
ied previously [l8l - |20| via both numerical and analytical 
approaches. Particularly when i/> 1, we have 



f(y)=log(2y)/(±y) + 0.13/y +0(y),. 



(10) 



The function K(y) can be determined numerically. To 
explore the functional form of K(y) and check the scaling 
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FIG. 3. Plots of function K(y) and f(y) based on compu- 
tational simulations. The dashed line corresponds to the ex- 
pression f(y) — a(y) found in Eg. (|13p that fits K(y) well. 



property proposed in Eq.®, we first run simulations for 
networks with different size L\ and Li but same shortcut 
range R and randomness p. The mean path length is £\ 
and £2, respectively. From on Eq.®, K can be found as 
(£% — I (L\ — Lq). The above process is repeated with 
various choices of L, R and p and each time we have a 
value of K. We then plot the K values found numerically 
as the function of y — rLp. If the scaling form holds, the 
plot of K(y) will collapse to a single curve, otherwise, the 
points will be scattered. 

In Fig [3], we show the K(y) curve obtained as dis- 
cussed above. As it displays a single curve as function of 
y, the scaling form in Eq. (|9j) is confirmed. To step fur- 
ther, we compared the function K{y) with f(y) and found 
that K (y) is very close to f(y) when y <C 1. When y>l, 
K(y) can be well fitted by the function K(y) = 2^+^l, 
which coincides with the second and other higher order 
terms in Eq. fTU]) . Based on these observations, we pro- 
pose the expression of I for large y values as: 
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(ii) 



which fits the simulation results quite well as shown in 
Figg]. 

In Eq. fTTl) . I changes logarithmically with L, which 
confirms our previous observation of small-world effect 
for large L values. More significantly, it indicates that 
the small- world effect appears when y = rLp 3> 1. It 
further indicates that global shortcut range is not nec- 
essary for small-world properties: even if the shortcut 
covers only a small fraction of the network, the small- 
world phenomenon will still occur as long as the network 
size is large enough. 

Eq. (JTTJ) gives the expression of £ in one limiting case. 
To find the full expression of £, we need the numerical fit 
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FIG. 4. Plots of numerical results of £ and Eq.fTT]) for 
L = 10,000,000, p = 0.004 and various values for r > 0.01. 
The choice of the L,p,r values guarantees the satisfaction of 
5 > 1. The blocks are the numerical results and the line 
corresponds to Eq.fTTJ. In the inset, we show the how term 
L(l — r)a(y) varies with network size L given r — 0.08 and 
p = 0.004. 



of K (y) for all y values. Based on the form of K(y) we 
found in the two limiting case, we related K (y) and f(y) 
in a way that 



K(y) = f(y)-a(y),. 



(12) 



a{y) is a positive function that approaches to when 
y — > and goes to log(2y)/(Ay) when y > 1. Using 
numerical fitting, we propose one possible form of a(y) 
as 



a(y) 



log(l 



24+y 



4// 



(13) 



which is displayed in Fig.[3]- 

The I in Eq.® can then be expressed as a combination 
of f(y) and a(y) as: 



£ = Lf(y)-L(l-r)a(y). 



(14) 



The first term is equivalently the mean path length of a 
SWN with rescaled randomness p' — rp, which always 
increases with L. The second term will be close to 
when L is small and then increase rapidly when L is 
large enough, as shown in the inset of Fig.(|4]). This ex- 
plains the nonmonotonic behavior of £ in SRSWN shown 
in Fig [2]. When L is small, the first term will dominate 
and £ will increase with L. After L reaches a critical 
value, the term La(y) will become non-negligible which 
can take the value of 1 down. This corresponds to the 
part where £ drops as L increases. Finally, when L is 
large such that rLp ^> 1, lb follows Eq. (fl~Tj) and in- 
crease logarithmically with L, similar to that of SWN. 
The difference of J between SRSWN and SWN derives 



5 



from the term a64( - 1 l^ 1 " 9 ^ , which yields large values 

for small r and p. This is why £ of SRSWN is larger 
than that of SWN even when L is large. (Note that the 
term °- 64 ( 1 4 p r + '° a ^^ can be negative for r > 0.5, which 
is mostly due to the error in numerical fit. But its ab- 
solute value is very small and we consider the error as 
acceptable.) 

III. SUMMARY 

In conclusion, we have studied a new variant of SWN 
model called SRSWN, where the shortcut length is con- 
strained within rN . The SRSWN shows the same degree 
distribution (see Appendix) and similar scaling proper- 
ties as SWN. However, it also exhibits distinctive struc- 
tural features which are not seen in SWNs. The most 
striking of these is the nonmonotonic behaviour of the 
avergae path length: for a range of system sizes, as 
the number of nodes increases the mean path length de- 
creases, i.e. the network appears to shrink as the number 
of nodes increases. It would be of interest to sec if this 
behavior is also reflected in more realistic networks which 
are expected to have a finite range of shortcut links. Our 
results also demonstrate that a global range for the short- 
cuts is not necessary to obtain the small- world effect. For 
SRSWNs, we derive analytical expressions for quantities 
of interest such as the mean path length which are in ex- 
cellent agreement with simulations. The results obtained 
thus suggest that similar approaches could be used to ac- 



curately characterize structural properties of diverse net- 
works which build on the original SWN model to include 
more realistic generalizations. 

The authors acknowledge funding support from NSF 
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IV. APPENDIX 

The distribution p(l) is defined as the probability that 
a node has I shortcuts attached to it. In SRSWN, p(l) 
amounts to 

P(0 = p(17)(i-^r™<^)<-' 

+ < i -rt( rl r 1 ) (i -^ ) ^ ,( ^ i5) 

The first term denotes the case such that there is one 
outgoing shortcut and I — 1 incoming shortcuts, while the 
second term is the case that all I shortcuts are incoming 
ones. 

In this paper, we focus on the case when R = rL value 
large. In this case, Eg. ([TBI) converges to a combination 
of two Possion distributions as 

p(l)= P j^e-v + (l-p)^e-P. (16) 

which is the same as that of SWN. 
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